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Abstract 

In this paper we consider the most general least-order derivative theory of gravity 
in which not only curvature but also torsion is explicitly present in the Lagrangian, 
and where all independent fields have their own coupling constant: we will apply this 
theory to the case of ELKO fields, which is the acronym of the German Eigenspvnoren 
des LadungsKonjugationsOperators designating eigenspinors of the charge conjuga- 
tion operator, and thus they are a Majorana-like special type of spinors; and to the 
Dirac fields, the most general type of spinors. We shall see that because torsion has 
a coupling constant that is still undetermined, the ELKO and Dirac field equations 
are endowed with self-interactions whose coupling constant is undetermined: we dis- 
cuss different applications according to the value of the coupling constants and the 
different properties that consequently follow. We highlight that in this approach, 
the ELKO and Dirac field's self-interactions depend on the coupling constant as a 
parameter that may even make these non-linearities manifest at subatomic scales. 
Keywords: Torsion tensor, ELKO and Dirac fields 
PACS: 04.20. Gz, 04.20.-q, 04.50. Kd, ll.10.-z 

1 Introduction 

The debate in general relativity about the presence of torsion, although one of the less 
known, is nevertheless one of the most fascinating, and several reasons have been put on 
the table, both against and in favour of torsion, ranging from philosophical to mathe- 
matical and phenomenological aspects of science, covering the foundational, theoretical 
and applicative domains of physics: the very first of these reasons is related to the fact 
that once the covariant derivative is written in the most general way, it is defined upon a 
connection that in general is not symmetric in the two lower indices, and therefore torsion 
does not generally vanish; the inclusion of torsion beside the curvature makes it possible 
to introduce a torsion-spin coupling beside the curvature-energy coupling, fully realizing 
the geometry-matter coupling prescription, which is lame otherwise; if we focus either on 
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the standard or on more exotic spinor fields, the torsional effects manifest themselves as 
spinorial self-interactions. Nevertheless, a problem that still needs to be addressed is that 
torsional effects appear to affect the dynamics only at the Planck scale, which would make 
torsion nearly negligible in almost every practical situation; this circumstance is due to the 
fact that, since both torsion and curvature have the same geometrical origin, then they are 
supposed to have the very same coupling constant, and a torsion having the gravitational 
coupling constant is what makes torsion interact weakly except at the Planck scale. 

On the other hand, despite torsion and curvature come in fact from the geometrical 
background, nevertheless they are independent fields, and therefore they should possess 
two independent coupling constants; previously [1] we have proposed a least-order deriva- 
tive generalization of the Hilbert-Einstein-Cartan-Sciama-Kibble HECSK Lagrangian in 
which a quadratic torsional term is present, with its own coupling. We have eventually 
studied an application to the case of the Dirac field. We have shown that when properly 
tuned there are consequences that could be detectable at subatomic scales. 

However, since torsion has three irreducible decompositions, then each of these should 
be thought as an independent field, entering with its own coupling constant; in this paper, 
we shall propound the most general least-order derivative generalization of the HECSK 
Lagrangian in which one quadratic torsional term for each of the three irreducible parts 
of torsion is present, and these three quadratic torsional terms will have three coupling 
constants. We will study applications to two cases: the case of ELKO spinors, whose name 
means Eigenspinoren des LadungsKonjugationsOperators, which designates eigenspinors of 
the charge conjugation operator, so that they are a Majorana-like spinor for which the 
problem of the mass term is solved by employing second-order derivative scalar-like field 
equations; the case of Dirac spinors, the most general type of spinor fields with least-order 
derivative field equations. In both cases the choice of finely-tuned couplings will have 
consequences stretching from cosmology, where we will study the case of ELKO fields, to 
particle physics, in which we will discuss the case of Dirac fields. 

The paper is organized as in the following: in the first section, we shall discuss why we 
propose a generalized gravitational Lagrangian for curvature and torsion, all independent 
fields having its own coupling constant, writing down the most general Lagrangian and 
the corresponding field equations; in the second section we shall apply this model to both 
ELKO and Dirac fields, showing how torsion endows ELKO and Dirac field equations with 
self-interactions having free coupling constants. 

2 The geometrical structure of the field equations 

In this paper, we consider spacetime (l+3)-dimensional manifolds possessing curvature and 
torsion, with its three irreducible decompositions, accounting for a total of four indepen- 
dent geometrical fields, given by the metric tensor, the torsion completely antisymmetric 
dual of an axial vector, the torsion trace vector and the torsion non-completely antisym- 
metric irreducible tensor; the metric tensor is given by g aa with the covariant derivative 
denoted by defined through the dynamical connection T^ a : metricity Dg = holds. 

In terms of the connection alone it is also possible to define two additional structures 
that are fundamental when studying the non-commutative properties of the covariant 
derivatives, one of which is the Riemann curvature tensor defined as usual by 

n p — b r p — 8 v p + v p r CT — r p r a (2 ~\ ) 

while the other is the Cartan torsion tensor defined as the antisymmetric part in the 
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connection as in the following 

Q f \v = ^ P fiu ~ ^ufi (2-2) 
so that once the metric is considered, also the contorsion tensor can be defined as 

and because of their definition together with the metricity condition Dg = 0, we have that 
the curvature G aripl , is antisymmetric in both the first and second couple of indices, while 
torsion Q ppv is antisymmetric in the last two indices and contorsion K pilv is antisymmetric 
in the first two indices; due to these antisymmetry properties, the curvature tensor has 
one independent contraction given by G p r)pv — G r) „ with contraction G r]v g TI1 ' — G as 
usual, while both torsion and contorsion have one independent contraction that is given 
by K vp p = Q p pv = V v : with trace vector V v and completely antisymmetric dual of the axial 
vector 2K pilv e ppva = Q P p V £ pllva — W a , it is possible to see that the most general torsion 
can actually be decomposed in irreducible decompositions according to the identities 

Q ppv — ^appu ~t~ 3 (Vi/Qpfi Vp9pv^) 

+ [Q P p V - \W a e appv - | {V u g P p - V^gpt,)] 

and also 

+ [K pup - YjW a e ap , vp - \ {V p g pv - V v g pp )] 
so that we may define the tensors 

^ppv Qppu ^appis 3 {Vu9pp V p g pv ) 

and also 

G p i/p — K pup ~Y}}^ E-apvp 3 (Vpfjpv VvQpii) 

as the non-completely antisymmetric traceless parts with which we have that 

Q ppu ^appu ~\~ 3 ^YvQpp Vp9pv) ~t~ Tppv 

and also 

K pvp — Y2^^ ^&pvp 3 (V p gpv VvQpfi) G pup (2.9) 

are the decompositions of torsion and contorsion in their three irreducible parts identically. 

Another notable geometrical attribute is the Levi-Civita covariant derivative V M de- 
fined via the Levi-Civita symmetric connection A£J CT ; its associated curvature is the Rie- 
mann metric curvature tensor R p vpl/ , with contraction R VL , contracted as R called Ricci 
metric curvature tensor and scalar as usual. The importance of the Levi-Civita torsionlcss 
connection is that with it, and the contorsion tensor, we have that the most general dy- 
namical connection can be decomposed according T p a = A£J CT + K^ aa which will be used to 
separate contorsion, or equivalently torsion, from all remaining purely metric quantities. 

Equivalcntly, we may employ the vierbein or tetrad formalism, in which orthonormal 
tetrad fields are introduced so that g a fs^a^b = r l ah -> where rj a b is the Minkowskian 
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(2.4) 

(2.5) 
(2.6) 
(2.7) 



(2.8) 



matrix, while the covariant derivative denoted with is defined through the dynamical 
spin-connection : conditions D£ — and Drj = are imposed, and when the former 
condition D£ = is explicitly written it reads T h m = Cj^p(^ p al/t + &^k) while the latter 
result into P 6 "',, = — r- 7 „ so that the first formula is used to establish a link between the 
connection and the spin-connection while the last antisymmetry represents a property of 
the spin-connection that will be important to describe the underlying Lorentz structures 
when we shall be dealing with Lorentz spinorial representations, in the following. 
In this tetrad formalism the curvature tensor is translated as 

= d a v^ - d«Ti a + r^rl - r^rl (2.10) 

while the Cartan torsion tensor is 

-Q\u = d^C - + r^i - r£4 (2.11) 

as it is easy to check: notice that one of the advantages of this formalism is to highlight 
the similarities between curvature and torsion as strengths of the spin-connection and 
tetrad- field potentials as a Poincare gauge theory in the Yang-Mills sense [5] . 

Before proceeding, we recall that it is possible to define a geometry of complex fields, 
where the gauge-covariant derivative is defined through the gauge-connection . 

From the gauge connection we can define 

F„ u = d^A v - d v A„ (2.12) 

as the strength of the gauge-connection potential as an internal U(l) gauge theory in the 
Yang-Mills sense, as it is usually done in any introduction of quantum field theory. 

Another of the advantages over the coordinate formalism of the tetrad formalism is 
that the transformation laws are no longer general coordinate transformations but Lorentz 
transformations, which have an explicit structure suitable to receive other representations 
such as the complex representation in which the inclusion of the complex fields above 
would actually fit perfectly: to see how, we consider first of all the introduction of the 
set of complex matrices j a of the Clifford complex algebra {7j,7j} = 217^ giving the set 
of complex matrices <J a b as [7^,7^] = 4<Tjj such that {~fi,Cjk} = i£ijkqH q and since in 
(1 + 3)-dimensional spacetimes all representations are equivalent, we only need to choose 
one and in the following we will always work in chiral representation; the set of er^ matrices 
is the set of complex generators of the infinitesimal Lorentz complex transformation, called 
spinorial transformation, for which the spinorial covariant derivative D p is defined by 
the spinorial connection the spinorial constancy of the 7^ and er^ matrices can be 
equivalently translated into the fact that the spinorial connection can be written as a 
Lorentz-valued connection plus an abelian term which can therefore be identified with the 
gauge-connection potential as T M = ^r ab ^a a b + iqA^I for any value of the q label. For a 
list of topological properties for such global constructions, see for instance references [3JH]. 

Now to see what these field equations actually are, we will employ the usual method 
of establishing on the basis of general principles a given Lagrangian and varying it with 
respect to the independent fields: the formalism through which a general Lagrangian 
depending on curvature and torsion as well as spin-connection and tetrad-fields has been 
established and intensively discussed in most of its consequences already in El El IE] 
and we are not going to discuss it here; the point we need to underline is that this 
formalism is the gravitational equivalent of what would be a formalism in which the general 
Lagrangian depends on the field strength as well as the gauge-connection potentials, as in 
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the usual Yang-Mills theory. Consequently, by establishing the proper matter Lagrangian, 
its variation with respect to matter fields will give the corresponding matter field equations. 

One of the advantages of both the method presented in references [2J 13] and the for- 
malism developed in references |H] is that according to these there is an equivalent 
role covered by both curvature and torsion in the underlying framework upon which the 
dynamics is built, and consequently a similar equivalent importance of curvature and tor- 
sion should be expected in the gravitational dynamics as well; therefore it should not be 
surprising that such a complementarity between curvature and torsion should be found in 
the construction of the gravitational Lagrangian of the theory. However this is not the 
case in general, because as in the simplest HECSK theory, the field equations come from 
the HE Lagrangian Cg = G clearly favouring curvature over torsion; consequently within 
the field equations torsion enters with the same coupling constant of the curvature, that 
is the gravitational coupling constant, and this is the reason for which torsional effects arc 
negligible for all practical purposes. At this point a reasonable question we may want to 
ask is whether the HECSK theory should not be regarded as a theory that is too simple. 

Inasmuch as both curvature and torsion are taken as independent geometrical fields, a 
more general treatment allowing torsion not only implicitly within the curvature but also 
explicitly on its own would appear to be more elegant: an example of theories in which 
torsion must be considered explicitly and not only implicitly through the curvature, so 
that the Lagrangian is constructed on an intertwined mix between torsion and curvature 
tensors, are the conformal theories of gravity [9]; however, these theories have higher-order 
derivative Lagrangians, and consequently their torsion propagates in regions in which the 
effects of torsion are expected to be vanishingly small [TO], and if torsion does actually 
vanish, so that the field equations of Weyl gravity are recovered [TJJ , there are nevertheless 
solutions that do not reduce to those of the known Einstein limit [12] H3j and therefore it 
is not surprising that discrepancies with observations do indeed arise |14j . This behaviour 
is not a peculiarity of conformal gravity, but of all higher-order differential theories of 
gravitation as a whole; a wise choice might consequently be to look for theories that have 
least-order derivative Lagrangians, so to recover the proper approximation. And among 
all of these theories, it is widely known that the HECSK theory is the simplest we have. 

So we have to face a delicate situation, in which we have to carefully find a balance 
between the push-forward, toward a generalization in which torsion is included beside cur- 
vature into the Lagrangian, and a pull-back, for which least-order derivative Lagrangians 
alone are considered in order to ensure the existence of the non-propagation and therefore 
the expected limit. Thus we may now ask: is there a more general HECSK-like theory we 
can consider? An answer has been recently given in [JJ , leaving however the possibility for 
another question: what is the most general HECSK-like theory we can possibly have? 

In order to look for this answer, we have to notice that a Lagrangian that is least-order 
derivative can only be linear in the curvature and quadratic in torsion: as an easy compu- 
tation shows, for a linear expression of curvature we can construct only one independent 
scalar, the Ricci scalar G, and for a quadratic expression of torsion we can construct three 
independent scalars, given by the Q v vp Q p pp , Qpp.vQ ppv , QpiivQ vp,p '- 1 the decomposition of 
the curvature in terms of the metric curvature plus torsion and the subsequent decompo- 
sition of torsion into its three irreducible parts shows that we may equivalently construct 
the Ricci metric curvature scalar R, with V^V^, W P W P , T ppu T ptJjU proving that all irre- 
ducible fields have been accounted. Therefore, a Lagrangian that will take into account 
the curvature scalar and the three torsional terms will correspondingly have to be settled 
in terms of four coupling constants, that is the Newton gravitational coupling constant 
and the three torsional coupling constants given by A, B, C independent on one another 
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and yet to be determined; such a Lagrangian will therefore be written in the form 

C G = G + AQ v V)1 Qw + BQ PI1V QP^ + CQ ptlu Q v ^ (2.13) 

as the most general that is linear in the curvature and quadratic in torsion, that is the most 
general least-order derivative Lagrangian defining the most general HECSK-like theory. 

This gravitational Lagrangian is supplemented by the electrodynamic Lagrangian given 
in terms of the charge q as the coupling constant of the Maxwell theory; such a Lagrangian 

C E = -kF a pF° P ( 2 - 14 ) 

is the only least-order derivative Lagrangian that we could possibly define. 

Finally the matter Lagrangian Cm in terms of the mass m of the matter field and with 
Planck constant h normalized to the unity will complete the Lagrangian of the theory. 

This complete Lagrangian, upon variation with respect to all independent field in- 
volved, will yield the field equations 

C (Q^ p - Q v P p + 2Qpp u ) + B (2Q y PP - 2Qp v p) + A {V v g p P - Vg p,/ ) + 

+ (Qpp v + VPgP" - V v g p P) = -S p P u (2.15) 

for the torsion-spin coupling and 

C (d^Qpp* - D^QPP a + V^Qpp* - V^Qpp" + Q Saa Qj - \Q e ™Q^eg pa ) + 

+B (2D^Q a PP + 2V tl Q a PP + 2Q e ° a Q e<T p - QP e<J Q a 0a - \Q 6 ™Qo a *g pa ) + 

+A {-D a V p + D^VPgP a + \V^VPg pa ) + 
+ (G pa - \Gg pa ) - \ {\g pa F 2 - FP e F a g ) = \TP a (2.16) 

for the curvature-energy coupling, together with 

\F^Qpp» + FPPQ^ + D a F°P = qJP (2.17) 

for the gauge-current coupling; the spin and energy densities S P,1V and T pa and the current 
density J p will be given along with the matter field equations when the Lagrangian of 
matter will be set explicitly: the Bianchi identities are converted by these field equations 
into the conservation laws given by the following expressions 

DpSPP" + V p Spp» + = (2.18) 

and 

DpTPP + V^Tpp - T^Q°PP + Sp^G^P + qJ F^ = (2.19) 

and also 

D P JP + V P JP = (2.20) 

as conservation laws for the spin and energy and for the current density, which will turn 
to be valid once the matter field equations will be implemented within the model. 

Before going forward, we must spend some words to study the general structure of 
the field equation for the torsion-spin coupling (|2.15p for which, as torsion or contorsion 
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can be decomposed in three irreducible parts, the field equations can correspondingly be 
decomposed in three irreducible parts: then we have one field equation for the trace vector 

aV v = ~SP pu (2.21) 

one field equation for the completely antisymmetric irreducible tensor 

bW a = -SP^e pilva (2.22) 

and one field equation for the non-completely antisymmetric irreducible tensor 

cT p^ = _ ^ S pv.v _ i s ^r £ ^ Ta£ap pu _ i (S a av g pli - S^gP")] (2.23) 

and after renaming 3 A + 2B + C — 2 = a, AC — 4B + 1 = 6, IB + C + 1 = c as new coupling 
constants for the sake of simplicity, we have several possible scenarios: 

1. The most general torsional case: in the case a^O, b ^ 0, c^O the model is the 

most general HECSK gravitation capable of accounting for all the parts of the spin of 
the matter field; this theory would be able to describe the most general higher-spin 
boson or fermion with no constraint whatsoever. 

2. Double-part decomposition of torsion: in either of the cases a = 0, b ^ 0, c^O 

or b = 0, a ^ 0, c ^ or c = 0, a ^ 0, b ^ the model reduces to a simpler 
HECSK gravitation without the possibility to couple, respectively, to the trace or 
the completely antisymmetric or the non-completely antisymmetric traceless part of 
the spin of the matter field; this theory would be able to describe more general higher- 
spin bosons or fermions although constraints over the spin are found also in this case, 
except for special instances in which particular symmetries are implemented, as for 
the Friedmann-Lemaitre-Robertson- Walker spacetime, where only the two vectorial 
parts of torsion are compatible with cosmological constraints, so that the case given 
by the choice a = 0, b ^ 0, c ^ 0, that is the first instance, fits perfectly. 

3. Single-part decomposition of torsion: in either of the cases b = c = 0, a ^ or 

a = c = 0, & ^0ora=i>=0, c^0 the model reduces to a simpler HECSK 
gravitation coupling, respectively, to the trace or the completely antisymmetric or 
the non-completely antisymmetric irreducible part of the spin of a given field; this 
theory would be able to include also more vector fields with mass and higher-spin 
tensor fields and spinor fields can be described, but apart from the second instance, 
a = c = 0, b 0, coupling to the simplest spinor field, the Dirac fermion, in more 
general, higher-spin bosons or fermions, constraints over the spin are found. 

4. The simplest torsionless case: in the case a = b = c = the model reduces to the 

simplest HECSK gravitation in which there is no coupling for any decomposition of 
the spin; this theory describes only fields whose spin density tensor S p,lu is equal to 
zero, and therefore it is able to account for scalars or vector fields coming from a 
gauge invariance and nothing more, so that more general vector fields with mass or 
higher-spin tensor fields or even spinorial fields are excluded from this scheme. 

For a thorough discussion about the decompositions of torsion see |TS] . Having discussed 
how torsion should be decomposed, it is easy to see that in all cases 4, 3, 2 the decom- 
position of torsion that corresponds to the vanishing coefficient cannot be determined 
unambiguously, and it is only in the case 1 that it is possible to unambiguously determine 
all parts of torsion and invert the field equations for the torsion-spin coupling (|2.15p in 
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order to get the expression of torsion in terms of the spin; then in the remaining field equa- 
tions (|2.16p and (j2.17[) it is possible to decompose all covariant derivatives and curvatures 
into the Levi-Civita purely metric covariant derivatives and curvature plus contorsional 
contributions, which can be written in terms of torsion and eventually substituted through 
the field equations above in terms of the spin density of the matter fields. 

In the case in which the spin density vanishes, then torsion and contorsion vanish as 
well, the gravitational and electrodynamic field equations are the HE and Maxwell field 
equations so that the Newtonian limit is recovered as desired. 

3 The coupling to fermionic fields 

In this section we apply the presented theory to a couple of fermionic fields given by spinors 
with spin-i; the special case of a Majorana-like spinor with higher-order derivative field 
equations, the ELKO field; and the more general spinor but with least-order derivative 
field equations, the usual Dirac field. For them the spin and energy densities together with 
the current will be given explicitly, and the matter field equations will be considered. 

3.1 The ELKO field 

A few years ago, a new form of matter field called ELKO matter field has been introduced 
in the panorama of physics: this form of matter field is represented by spinors with spin-i 
A verifying the charge conjugation conditions 7 2 A* = ±A for self- and antiself-conjugated 
fields respectively, and thus they are a type of Majorana matter field [TB]; with this defini- 
tion, ELKOs are clearly topologically neutral fermions, and because a topological charge 
is what keeps localized the otherwise extended field, it follows that the absence of such 
charges ensures that nothing protects the field from spreading, allowing them to display 
non-local properties: as a consequence, according to the Wigner prescription, for which 
fundamental fields are classified in terms of irreducible representations of the Poincare 
group, ELKOs belong to a non-standard Wigner class [17] . Thus defined, one has that for 
the ELKO matter field to possess mass, without the introduction of Grassmann-valued 
fields, they have to obey second-order derivative matter field equations, and therefore they 
must have mass dimension 1 [T5J [TO] ; however, still referring to the Wigner classification 
scheme, fundamental fields are labelled in terms of both mass and spin quantum numbers, 
implying that matter fields have both energy and spin density tensors, and thus requiring 
both curvature and torsion, if ELKOs have to include the most general coupling to gravita- 
tion [20, 21 . All of the properties found in the ELKO matter field definition contribute to 
some effect in their dynamics: for instance, the ELKO matter field is proven to have q = 0, 
that is it has no coupling to electrodynamics and so it is a natural candidate for dark mat- 
ter [22] , while their completeness relationships indicate that there must be a preferred axis 
of locality, probably having something to do with the privileged direction that arises within 
the scheme of Very Special Relativity [23] ; applications to cosmology, in particular dark 
energy, have been addressed in different works [23] [23 [23 [23 12H1 1211 130] ; their algebraic 
structures have been discussed quite thoroughly [32 1321 1331 133] , whereas their dynamics 
in terms of exotic properties has been studied as well [33 [33 GUI I3H] ■ An alternative 
approach starting from quantum field theory formalism has been pursued [391 140] , while 
a variety of phenomenological applications have been investigated [411 EE] S3 EH S3 [46] . 

Now, that ELKOs have second-order differential field equations and a coupling to 
torsion has two important consequences: first, the ELKO field equations have kinetic 
term with two derivatives, and secondly, their spin density tensor is differentially related 
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to the torsion tensor; these two consequences taken together imply that, within the ELKO 
field equations, there is the appearance of derivatives of torsion, which is itself containing 
derivatives of the ELKOs, and thus additional second-order derivatives of the ELKOs 
do arise, which may induce acausal propagation and in the high-energy ultraviolet limit 
singularities may form. Nevertheless, all torsional fermionic back-reactions are proven to 
cancel exactly, ensuring the ELKO causal propagation while in the high-energy limit the 
ELKO gravitational pull tends to vanish, showing that ELKO matter fields are causal and 
they have a gravitational asymptotic freedom for which their topological non-locality is 
extended to include a dynamical form of gravitational non-locality [UJ HH1 US] • ELKO 
matter has then been generalized, both in terms of the ELKO matter field dynamics itself 
and in terms of the coupling to the gravitational background in various ways [501 1511 152] . 
The matter fields defined as ELKO and ELKO dual are explicitly given as 

7 2 Ai T = 7?A ±T X* ±T 7 2 = -V X±t (3-1) 

with relationship between ELKO and ELKO as A=p±= ±iAj_ =(= 7 in which 77 = ±1 for self- 
or antiself-conjugate fields and with the label ±=p indicating that the fields decompose 
into irreducible chiral projections that are eigenstates with positive/negative or nega- 
tive/positive eigenvalues of the helicity operator: ELKO and ELKO dual have explicit 
decomposition in terms of the irreducible chiral projections given by the following 

A ±T = ( ±z^,x 2 T iL% ) A ±T = ( _ v L J L , ± ) (3.2) 

with the label ± designating the eigenstate with positive or negative eigenvalue of the 
helicity operator respectively given by 



X+-= ( ~m -iq*) A+_ 



A- + = (0 - m iq* ) A_ + 



when their spin and momentum have the same direction; all ELKO special algebraic 
properties come from their special structure in terms of Majorana-like spinors, while there 
is nothing special in their differential properties as these depend only on the fact that they 
are spin-i spinor fields. Given the specific algebraic structure and the usual differential 
structure we may now proceed and write the most general ELKO Lagrangian. 

The Lagrangian we will employ will be the most general we can possibly write, and 
since ELKO fields are spinors of mass dimension 1 this consequently means, first of all, that 
the Lagrangian for ELKO must account for all possible terms built with two derivatives, 
and finally, that for ELKO there are two such terms: the first is the one we would have 

have in the simplest case D a A D a X but also the term D a \ cr a ^ DpX can be taken since 
due to the presence of torsion it does not reduce to the divergence of a vector, negligible 
in the action; therefore, the most general of such Lagrangians for ELKO will eventually 
be given with both these terms, and one normalization factor aside, it will contain one 
parameter p to be determined. The most complete Lagrangian for the ELKO field is then 



/ d \ 



d 



irjd* J 
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d 

irjd 
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(3.3) 



(3.4) 
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given by the following form 

£elko = D a A (v a0 I + p<r afi )Dp\ + m 2 A A (3.5) 

in terms of the mass m of the matter field. 

Its variation gives the spin and energy densities 

Sfiap — \ (yDfi A cr a p\~ A (TapDf+Xj + 

+| (D p A (T p n<TapX- A <T a p<T^ p D p Xj (3.6) 

and 

T„ v = (Dp A D u \ + £>„ A D^X - g„ v D p A D p a) + 
+p (.D„ A (r w D p \ + D p A o> M D„A - g^Dp A o-^D^a) + g^m 2 A A (3.7) 

without field equation for the current density since ELKO fields are neutral and while the 
matter field equations are 

(£> 2 A + V^D^X) + p (a-P p D p D^X + VpCr'^D^X) + m 2 A = (3.8) 

to be taken into account: as soon as the matter field equations are employed the energy 
and spin densities undergo the conservation laws (|2.18p and (|2.19[) above. 

The whole system of field equations is thus given by the field equation for the torsion 

-2C {Q pvp - Q vpp + 2Q« ll/ ) - 2B {2Q vp P - 2Q p - vp ) - 2A (V v g pp: - V p g pv ) + 
-2 (QP pu + V»gP v - V u g pp ) = (d p A er^X- X ct^D p \^ + 

+p (p^ X cr^Pa^X- X a^a^D^Xj (3.9) 

and the curvature 

2C (DuQw - D^Qpp* + V^QPP a - V^QPP a + Q 6aa Qj - \Q 6a * Q^eg pa ) + 
+2B (2D^Q a PP + 2V t j i Q a PP + 2Q e ° a Q ea p - Q p6 °Q a Ba - \Q ean Qe^9 pa ) + 
+2A (-D a V + D^V^gP* + ^V f ,V p g pa ) + 
+ (2G pa - Gg pa ) = (p p A D a X + D a X D P X - g pa D 7r X D n X) + 

+p (D a X tr^D^X + D n X (T*PD a X - gP a D n A cr™D a X) + gP a m 2 X X (3.10) 

along with the matter field equations 

(D 2 A + FD M A) + p (o-PPDpD^X + V p a p ^D^X) + m 2 X = (3.11) 

which are rather complicated, and the decomposition of the torsion-spin coupling field 
equations, let alone the decomposition of torsion and its substitution with the spin in all 
other field equations, may be very difficult to perform. 

However, we may employ the freedom this theory gives us to choose special fine-tunings 
for all parameters involved, which considerably simplify the solution of this problem, and 
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the following treatment; this is due to the fact that in such cases, in which some of the 
parameters a, b, c vanish, the corresponding trace, completely antisymmetric and non- 
completely antisymmetric traceless decompositions of torsion are not uniquely defined, 
and they may therefore be chosen to vanish without any loss of generality. In the next 
sections we are going to show how these simplification will be performed in two special 
examples in cosmological applications of particular interest. 



1. The FLRW-compatible ELKO field with % A 7 A = constraint. To study 
cosmological applications, we have to implement the cosmological principle, for which 
the metric is Friedmann-Lemaitre-Robertson- Walker metric, which will be chosen to be 
spatially flat as 



9tt = 1, 



9jj 



x,y,z 



(3.12) 



where torsion has only the time component of the vector parts V t and W t ; the ELKO field 
that turns out to be compatible with such symmetries may be chosen to be given by 



A+_= V2ip ( r] i ) A+_ = V2ip 



A-+= \/2> (Or) - i ) A_ + = y/2<p 



( 1 \ 





V -iri J 

( °\ 

1 

if] 

V oy 



(3.13) 



(3.14) 



with the additional property for which i X 7A = identically: in this way torsion has only 
the time component of the vector part Vt if we assume p = — 1 for the ELKO matter. 

In this case it is possible to see that the field equations for the torsion-spin coupling 
are decomposed to the single independent field equation for the torsion-spin trace vector 

aV v = d v (3.15) 

so that the field equations for the curvature-energy coupling are 

-2 (2£a) (D a VP - D^VfgP a - ±V^VPgP a ) + 
+ (2GP a - GgP") = [d p A D a \ + D a A D P X - gP a D 7r A D^X^j + 
- (p a A crP n D n X + -Djr A cr n PD a X - g pa D 1T X cr 7rcT D a X^ + gP a m 2 X 2 (3.16) 

and we have now to separate torsion away everywhere; after this is done, the field equations 
for the Ricci tensor are 



2RP° = f (V a VP + \V^VPgP a + \VPV a - |V%£T*) + \VPV a X 2 + 
+\Vp (V p A a^ a X- X (r (Sa V p X+ A ct^ct^V^A - V T A a^Pa^X^ + 

+ (V p A V"A + V" A VA - V" A tr^V^A - V T A it'Ta) - g pa m 2 X 2 (3.17) 

in which now torsion has to be substituted in terms of the ELKO fields: and after this is 
done as well, the time-time component will be such that only one term will contain the 
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parameter a so that for a very small this term will be much larger than all others leaving 

f«-Hw) 2 (3-18) 

describing the evolution of the scale factor a as a function of the cosmological time. 

As it is easy to see for a that approaches zero from negative values we have that 
the acceleration a is positive and very high, therefore this may turn out to give relevant 
contributions for the accelerated expansion of the universe. 



2. The FLRW-compatible ELKO field with A A = constraint. We continue the 
study of cosmological application by still considering the Friedmann-Lemaitre- Robertson- 
Walker spatially flat metric 



5oo 



1, 



9jj 



J = x,y,z 



(3.19) 



where torsion has only the time component of the vector parts V t and W t ; however in this 
case the ELKO field will be chosen to be 



A+_= V2if ( if} 1 ) A+_ = V2if 



A-+= V2<p ( irj - 1 



( 1 \ 





V -ivj 

( °\ 
1 

ir\ 

V oy 



(3.20) 



(3.21) 



with the additional property for which A 2 = identically: in this way torsion has only the 
time component of the vector part W t if we assume p = — 1 for the ELKO matter. 

The field equations for the torsion-spin coupling are decomposed to the one independent 
field equation for the torsion-spin completely antisymmetric dual of the axial vector 
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(3.22) 



so that the field equations for the curvature-energy coupling are 

(i^) (D^Wge^P 01 + \W a WP + ±W 2 gP a ) + 
+ (2GP a - GgP") = [d p A D a \ + D a A D P X - g pa D n A D*a) + 

- (p a A o-p^D^X + A cr n PD a X - g pa D 7r A a™D a X} (3.23) 

and we have now to separate torsion away everywhere; after this is done, the field equations 
for the Ricci tensor are 

2RP a = | {\7^W g eP e P a + \W a WP - \W 2 g pa ) + ±W p V a (i A 7A) + 



+ rfeWp [V p A <T^ a jX- A jcr^ a V P X+ A Tcr^cr^V^A - A cr^cr^jX 



+ (V p A V Q A + V" A V"A — V" A o-'* V ff A - A cr^V a A) 



(3.24) 
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in which now torsion has to be substituted in terms of the ELKO fields: and after this is 
done as well, the time-time component will be such that only one term will contain the 
parameter b so that for b very small this term will be much larger than all others leaving 

f«-^(w) 2 (3-25) 

describing the evolution of the scale factor a as a function of the cosmological time. 

And for b that approaches zero from negative values the acceleration a is positive and 
very high, turning out to be relevant for the accelerated expansion of the universe. 

It is intriguing that in these two complementary examples we eventually get an anal- 
ogous result at the end: of the four contributions to the ELKO energy, whether they are 
the standard or the additional one proportional to the parameter p, whether in terms of 
the torsionless derivative and the torsional potentials, the one coming from the additional 
energy proportional to the parameter p in the torsional potentials, that is the one that has 
the most genuine torsional origin, is the only one that remains to be relevant when the 
remaining parameter is chosen to approach zero from negative values; once this is sorted 
out, it controls the behaviour of the accelerated expansion of the universe. 

3.2 The Dirac field 

One of the reasons for which torsion was long neglected is to be found in the historical 
misunderstanding that the implementation of the principles of equivalence and causality 
was supposed to force the symmetry of the connection and consequently implying the 
vanishing of the torsion itself [531 [54] ; this misunderstanding however may be resolved by 
a deeper analysis, and if care is taken it is straightforward to acknowledge that even in 
the case in which the principle of equivalence holds then torsion is only restricted to be 
completely antisymmetric [55, SSI EZ1 EH1 ES] ; it is therefore quite intriguing that if torsion 
is completely antisymmetric then the spin must be completely antisymmetric as well, and 
the only spin that is completely antisymmetric without any constraint imposed is that of 
the simplest spinor field [501 US] ; this simplest spinor field is indeed the Dirac field, for 
which the completely antisymmetric dual of an axial current appearing as a self-interaction 
in the Dirac field equation is able to give a dynamical explanation of the exclusion principle 
as stated by Pauli [62j . The study of the Dirac equation is for this reason quite interesting. 
The Dirac field Lagrangian is notoriously the only least-order derivative Lagrangian 

ADirac = % -{^D^ - D^Y^P) ~ (3-26) 

in terms of the mass m of the matter field. 

Its variation gives the completely antisymmetric spin and energy densities 

S fia p = ^{j fi ,(T a p}ip (3.27) 

and 

T m = | (4>l^D a il> - Djfoytf) (3.28) 

and the current 

J, = (3-29) 

as soon as the matter field equations 

i^D^ip + '-YV^ - = (3.30) 
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are taken into account, and where the matter field equations also imply that the energy 
and spin densities and the current obey the conservation laws f|2 . 18[) and (|2.19l) above. 

The whole system of field equations is thus given by the field equation for the completely 
antisymmetric torsion 

bW a = \^Yl^ (3-31) 

and the curvature 

(tt) {D^W a e^ pa + \W a W p + j^W 2 g ap ) + (G pa - \Gg pa ) - 

-\ (}g pa F 2 - F pe F a e ) = i(<tpyD a ip - D a ^YiP) (3.32) 

and also the gauge field 

±F^W a e p ^ p + D a F ap = q^Yij (3.33) 

along with the matter field equations 

iY D n4> - mip = (3.34) 

for which the conservation laws (|2.18[) and (|2.19l) above will turn out to be important for 
the decomposition in torsionless and torsional quantities within all the field equations. 

After such a separation and upon substitution with the spin density of the Dirac field, 
we get the field equations for the Ricci tensor 

Rpa ~ 1 (\9 pap2 ~ F p9 F a g ) = § (^(''VV - V^ipYH) ~ \g pa miP 2 (3.35) 
and the gauge field 

V a F ap = #7^ (3.36) 

along with the matter field equations 

iYV^ + ^YlH^ -mip = (3.37) 

showing that the gravitational field equations decomposed into their symmetric part writ- 
ten in terms of the Ricci tensor and the electrodynamic field equations are formally those 
we would have in the torsionless case while the matter field equations become formally 
identical to what we would have had in the torsionless case with additional potentials 
of self-interactions having an independent coupling constant, and where the gravitational 
field equations in their antisymmetric part are implied by the spin conservation law (|2.18D . 

Focusing on the non-gravitational field equations, we see that they may be written 
after a Fierz rearrangement as 

V a F ap = qiPYiP (3.38) 

together with 

iY - •m^H'H^ — mip = (3.39) 

showing that the vector ipY" , P accounts for both electrodynamic and self-interaction. 

A further simplification may be given in terms of another Fierz rearrangement while 
having the coupling constant renamed as A = to get 

V M -0 - A (■0V - "07^7) - m i> = (3.40) 
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as the matter field equations will now study in specific situations. 

We remark that in the present theory as in the HECSK theory, the Dirac field equations 
are of the Nambu-Jona-Lasinio type [52]; however, in this model the Dirac field has self- 
interactions whose coupling constant is not determined [53] ■ It is important to notice that 
no matter what this value is, it would at first be thinkable to normalize it to the usual 
value of the Newton constant through the normalization of the Dirac field, but a deeper 
analysis shows that this process would also have the effect of changing the scale of the 
energy density within the gravitational field equations and it is thus unacceptable. 

So far as our knowledge is concerned, the theory we have here constructed is the most 
general in which the HECSK theory for the Dirac field can be generalized in order for the 
system of field equations, once decomposed in terms of torsionless quantities plus torsional 
contributions written as Dirac field self-interactions, to provide a free coupling constant 
for the self-interactions of matter while leaving the gravitation unmodified. 



3.2.1 The Dirac field self-coupling and its applications: 
Leptons with weak interactions and Higgs field 

To begin with, we take into account the case of a coupled system of Dirac spinorial matter 
fields, one of which is a spinor with mass while the other is a semispinor and therefore only 
compatible with the massless configuration; these spinors will possess spin-torsion coupling 
giving interactions between the two fundamental fields while in the standard model the 
couple of electron and neutrino may only have weak interactions: the aim is therefore to 
compare this case with that of the weak forces and see whether or not similarities arise. 

In this case the field equations for the electron and neutrino fields e and v are given 
by the previous field equations, but because now the system is constituted by a couple of 
spinors then the spin will be the sum of the two spins and the spinorial field equations are 

«7 M V M e — A (e^f^e-y^e + F7 ^Yl e ) ~ me = (3-41) 
iYV^v - Xej^jeYiy = (3.42) 

in which the fact that the electron is massive while the neutrino is massless is the reason 
that prevents these two fields to mix into a doublet; equivalently we have 



ry^V^e + 2\(cos9) 2 e / ye~fe - 
To 

*'7"V^ + j^gZ^v + ^W^e L = (3.44) 



+gtan#Z M 7^e - ^JL-Z^cl + -%W*^v -He-me = (3.43) 



once we define 

Z» = -A [2 (sin 6) 2 eYe - clY^l + V^v] (j 2 ^-) (3.45) 

^ = -M^V) (3-46) 
H = Aee2(cos6») 2 (3.47) 

showing that field equations (|3. 431) and (|3. 441) are formally identical to the system of field 
equations for the lcpton fields after the symmetry breaking in the standard model, al- 
though both weak and Higgs boson fields are here composite [64 . This is an important 
discrepancy with respect to the standard model, because, in our approach, we expect both 
weak and Higgs bosons to display internal structure whenever the energy is high enough 
to probe their potential compositeness, while, in the standard model, they are supposed to 
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be structureless at any energy scale. Similar results are obtained for leptons in Higgsless 
models and therefore without symmetry breaking [65] ; for leptons and hadrons in circum- 
stances in which the symmetry was unbroken see references |66l I67| . The phenomenon of 
oscillations for neutrino in the massless configuration have been studied in reference |68j . 

In order to reproduce the weak interactions the tuning A = Gf ensures the coupling 
constant to have the value of the Fermi coupling constant while in the case of oscillations 
the constant A must be tuned to the oscillation length: the field equations that have to be 
solved in order to get the relationship between A and the oscillation length are non-linear 
and therefore numerical analysis must be used. But even though it is unlikely that the 
experimental value for A in this case may fit that of the Fermi model. 

3.2.2 The Dirac field self-coupling and its applications: 
the Dirac field modelling condensed state physics 

We will now write the Dirac field equation for the single field in the standard representa- 
tion, the one that allows the slow-speed weak-field limit, in which the spinor field maintains 
only the large component given by the semispinor <f>, and for which the spinor field equa- 
tion is split into two field equations of which one reduces to zero while the other is the 
field equation for the semispinor (f> given by 

+ ^<r k VkCT a V a - A (fitf,) - 7710 = (3.48) 

compatibly with the fact that the spinor actually has only two degrees of freedom, and 
representing cold-matter fields with self-interactions with the structure of a mass term. 

Thus in the present theory we have that the Dirac field equations have the form of the 
Nambu-Jona-Lasinio field equation further approximated to the Pauli-Schrodinger field 
equation of the Ginzburg-Landau type |62) . In it the coupling constant A is to be set on 
the value of the specific condensed state system we would eventually like to study [53] . 

From now on we will neglect the electrodynamic field, which in particular means that 
the magnetic field will be vanishing, thus permitting the decoupling between the spin 
degrees of freedom, so that the semispinor will correspondingly decouple in the two com- 
ponents given by the scalars u and v, and after the field equations have been separated as 
well, the u and v field will have to obey the coupled system of field equations 

i^ + ^V 2 u~\u 2 u-(m + \v 2 )u = (3.49) 
+ ^V 2 v - Xv 2 v - (m + Xu 2 )v = (3.50) 

representing cold-matter fields after a process of bosonization, with self-interactions and 
interactions between the two components; notice that in the approximation for which for 
one field the other field's square may be considered constant, then the residual interaction 
becomes a correction to the mass term, and the two equations may be treated separately. 

Thus in the present theory we have that the field equations are approximated to the 
Pauli-Schrodinger field equation of the Ginzburg-Landau type and then decoupled to the 
Schrodinger field equation of the Gross-Pitaevskii type. In it the coupling constant A has 
to be set according to the value of the specific condensed state system we would eventually 
like to study and as it is clear it may be attractive as well as repulsive as we shall see. 

If we now focus on the field u alone, it is possible to see that for stationary states the 
field equation reduces to 

^V 2 u - Xu 2 u + eu = (3.51) 
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with total energy e either positive or negative, and whose general solution may be given 
according to the sign of the constants and parameters involved, therefore we have a total 
of four different cases: for attractive potentials we will set mX = —g 2 while for repulsive 
potentials we will set m\ — g 2 ; and in each case, negative energy will be written in terms 
of 2me = —lu 2 and positive energy will be written in terms of 2me = uj 2 in what follows. 

Now an interesting application is that of the unidimensional potentials in the coordinate 
x, as in this case exact solutions are easy to find: for attractive potentials the energy may 
only be negative and a solution is given by 

«<*) = ! safe ( 3 - 52 ) 

whereas for repulsive potentials we still have the case of negative energy with solution 

«(*) = i Sh^fe ( 3 - 53 ) 

and that of positive energy with solution 

u(x) = (3.54) 

V / g sin (uix) v / 

as it can be checked by direct substitution; it is easy to see that while for attractive 

potentials we always have bound states that are finite and tend to zero at infinity, for 

repulsive potentials the bound states diverge as \ but still tend to zero at infinity whereas 

i .... 

the unbound states diverge as ^ and exhibit a periodic behaviour. If we focus on repulsive 
potentials we have that alternative solutions are obtained for negative energy as 

„ (!C ) = -^tan(^) (3.55) 

and for positive energy as 

u{x) = _|_ tanh (-) (3.56) 

again as it can be checked directly; in this case there is the tendency to spread at infinity, 
since what is supposed to be the bound states exhibit a periodic behaviour while the un- 
bound states display the solitonic behaviour. The repulsive potential with positive energy 
is interesting since a general solution can be found in terms of the Jacobi elliptic functions, 
and because torsion should have the behaviour of an angular momentum giving rise to 
potentials with the structure of centrifugal barriers and the energy is usually assumed to 
be positive, and therefore it represents the most reasonable physical situation. 

However, for a complete study of the unidimensional potential, suitable boundary 
conditions must be assumed, and therefore we will next turn our attention to the unidi- 
mensional potential box with width I: we focus on the case of repulsive potential with 
positive energy supposing that the repulsive potential is much smaller than the energy so 
that g is much smaller than w; an approximated solution is given by 

u n (x) = sin (oj n x) + ^2 [yg sin (3w n x) - |w„a;cos (w n x)] (3.57) 
with boundary condition u(0) = u(l) = giving discrete uj n and thus energy levels e„ as 

^ = n 2 ^ + \\ (3.58) 

which is therefore clearly quantized, as expected for any bound state with such a set 
of boundary conditions; notice that even in the ground-state configuration the torsional 
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interaction generates a zero-point energy that is not equal to zero. This is quite intriguing, 
and further calculations should be performed to get higher orders of approximation. 

Unidimensional plane waves may also be employed to evaluate the general form of the 
energy levels for the free system to be 

e = m + ^ + Xu 2 (3.59) 

with u 2 being the density of the field itself, which in the interpretation of a collective state 
of plane waves represents the number of waves per unit volume u 2 = ^; on the other hand, 
it is also possible to go a little further and employ a thermodynamic interpretation of the 
kinetic energy as the temperature T of the system: considering both interpretations, the 
energy levels are reinterpreted as 

e = T + Xy (3.60) 

which is the expression for the energy level of the Van der Waals gas. By employing the 
equation of the energy (Jy) T = (^f) v T — P it is easy to see that the equation of state 
gives the pressure P as a function of volume and temperature as 

P-X^~ T ^j = (3.61) 

as the Van der Waals equation of state with generalized volume factor F(V) and with 
coupling constant A: for negative values the Van der Waals pressure gives an attraction, 
for positive values the Van der Waals pressure gives a repulsion. In this reinterpretation 
the matter field is seen as a classic wave distribution behaving as a Van der Waals gas. 



4 Conclusion 

In this paper, we have constructed a generalization of the Einstein-Cartan gravity in which 
both curvature and torsion are considered each with its own coupling constant, the one 
related to the curvature being the gravitational constant while the three related to torsion 
are still undetermined; we have applied this geometrical background to the case of Dirac 
and ELKO fields showing that the field equations are endowed with self-interactions whose 
coupling constants are yet to be determined: we have discussed different applications in 
which these self-interactions play different roles according to the value of the constant 
we assign, whether they are the ELKO field used in the standard model of cosmology 
to fit the accelerated expansion of the universe or the Dirac field both in the standard 
model of particle physics to fit the behaviour of fermions with weak interactions or for the 
condensed state approximation. We have insofar witnessed that the presence of torsion 
with independent coupling constants may either give rise to new forms of interactions, 
in cosmology, or to interactions that are eventually shown to have the same structure of 
already known interactions and potentials which have always been ascribed to something 
else, whether the context was the domain of particle physics or condensed state physics. 

A problem we have already discussed in [T] was related to the fact that the approach 
followed there had only one additional coupling constant which might have been unable to 
have the whole spectrum of possible applications fit into a single framework; the approach 
followed here has three coupling constants, and so freer in its structure: for instance, it is 
possible to consider both ELKO and Dirac field in a single Lagrangian, in a FLRW universe 

and choosing a specific form for the ELKO so that i \ 7 A — we get that the torsion 
only has the trace part sourced by this ELKO field and the completely antisymmetric 
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dual of the axial vector sourced by the Dirac field. Then it is possible to assign the values 
to the coupling constants as to have c = in order to ensure the compatibility with the 
symmetries of the background, a approaching zero from negative values so that ELKO 
may fit the accelerated expansion of the universe and b approaching zero from positive 
values to have Dirac fields with repulsive self-interactions able to give rise to the weak 
forces or bosonization in condensed states. This a posteriori shows that it is reasonable 
to consider the three different parts of torsion independent as they are sourced by fields 
that are totally unrelated; in particular, this implies that it is reasonable to assign three 
different coupling constants, and employ them for different purposes. But this also shows 
that this freedom is anyway rather limited, as for a given field the constant is assigned 
once and for all, and it is unlikely that the very same field may have self-interactions 
with strengths that differ from one situation to another. Therefore the discussion outlined 
in our previous work is here improved for what regards the possibility to set different 
strengths for different parts of torsion, whether we are in cosmology or quantum physics, 
but not so much as to set different strengths for the same part of torsion in particle physics 
or condensed state physics, and for that an even more general theory of gravity may be 
able to provide a running coupling constant, possibly scaling with the energy, such that all 
applications above, and maybe more, can fit into a single scheme, as anticipated in [69] . 

If this approach really works, torsion would not only be observable, but its effects might 
have already been observed, although we have not been able to recognize them as coming 
from torsion; and more effects may be on the verge of being discovered. In writing this we 
are aware of the fact that torsion may not be the answer for all problems of physics, but 
it possibly is for some of them, and our theory can allow us to see how. 
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